Climate change did not alter the effects of Bt maize on soil Collembola in northeast China

Bt maize is being increasingly cultivated worldwide as the effects of climate change are increasing globally. Bt maize IE09S034 and its near-isogenic non-Bt maize Zong 31 were used to investigate whether climate change alters the effects of Bt maize on soil Collembola. Warming and drought conditions were simulated using open-top chambers (OTC), and their effects on soil Collembola were evaluated. We found that the maize type had no significant effect on Collembola; however, the abundance and diversity of Collembola were significantly higher in the OTC than outside at the seedling stage; they were significantly lower in the OTC at the heading and mature stages. The interactions of the maize type with the OTC had no effect on these parameters. Therefore, Bt maize had no significant effect on soil Collembola, and the effects of climate warming and drought on soil Collembola depended on the ambient climatic conditions. When the temperature was low, collembolan abundance and diversity were promoted by warming; however, when the temperature was high and the humidity was low, collembolan abundance and diversity were inhibited by warming and drought. The climate changes simulated by the OTC did not alter the effects of Bt maize on soil Collembola.

www.nature.com/scientificreports/ Bt maize under the predicted climate future conditions of warming and drought will be helpful for scientifically evaluating the potential ecological risk of Bt crops to soil ecosystems under future climate change conditions. In the present study, Bt maize IE09S034 and its near-isogenic non-Bt maize Zong 31 were used as the experimental crops, and open-top chambers (OTCs) were used to simulate warming and drought, which are expected with climate change. We analyzed the effects of the maize type, OTC treatment, and sampling stage on the abundance and diversity of Collembola. Our study aimed to clarify whether warming combined with drought will alter the effects of Bt maize cultivation on soil Collembola. We hypothesized that (1) Bt maize would have no significant effect on the soil collembolan community; (2) climate change simulated by OTC would have a significant effect on the soil collembolan community; and (3) climate change simulated by OTC would alter the effects of Bt maize on soil Collembola.

Results
Soil parameters. The mean soil temperature was approximately 25 °C at the seedling stage, increased to approximately 35 °C at the heading stage, and then decreased to approximately 18 °C at the mature stage. Oneway analysis of variance (ANOVA) showed that for the same maize type, the mean soil temperature was ~ 1 °C higher in OTC than outside for all the sampling stages (df = 1, 7; P < 0.05). The soil humidity and root biomass were low at the seedling stage and were not significantly different between the OTC and outside (df = 1, 7; P > 0.05). However, with time, the soil humidity and root biomass in OTC were both significantly lower than those observed outside at the heading and mature stages (df = 1, 7; P < 0.05) (Fig. 1).

Collembolan community composition.
A total of 19,581 individuals belonging to 7 families, 13 genera, and 17 species of Collembola were sampled. In total, 4561 individuals were in Bt maize in OTC, 4660 were in non-Bt maize Zong 31 in OTC, 5141 were in Bt maize outside OTC, and 5219 were in non-Bt maize Zong 31 outside OTC. The dominant species in all treatments was Hypogastrura yosii (52.84% in the Bt maize in OTC, 49.27% in the non-Bt maize in OTC, 43.18% in the Bt maize outside OTC, and 44.93% in the non-Bt maize outside OTC) and Thalassaphorura macrospinata (14.84% in the Bt maize in OTC, 18.37% in the non-Bt maize in OTC, 25.77% in the Bt maize outside OTC, and 26.00% in the non-Bt maize outside OTC). Other relatively abundant species were Folsomides parvulus, Isotomodes sp., and Entomobrya sp.4, with dominance more than 3% and less than 10% in each treatment. The remaining species were relatively rare ( Table 1).
Impacts of maize type, OTC, and sampling stage on Collembola. The three-way ANOVA results showed that neither the maize type nor the maize type with other factors (OTC and sampling stage) had a significant effect on the collembolan abundance and diversity, indicating that transgenic cry1Ie maize did not influence the soil collembolan community. The sampling stage had a significant effect on the collembolan abundance and diversity. The OTC had no significant effect on the collembolan diversity, but had a significant effect on the collembolan abundance, and the interactions of the OTC with the sampling stage had a significant effect on the collembolan abundance and diversity (Table 2). Further ANOVA results showed that for the same maize www.nature.com/scientificreports/ type, the collembolan abundance, species richness, Shannon-Wiener index, and Pielou's evenness index were all significantly higher in OTC than outside at the seedling stage; however, all these parameters were significantly lower in OTC than outside at the heading and mature stages (Fig. 2), indicating that the effect of OTC on Collembola was related to the sampling stage, when the temperature is low (seedling stage), Collembola are promoted by OTC, and when the temperature is high and the humidity is low (heading and mature stages), Collembola are inhibited by OTC.

Discussion
We investigated the impact of the maize type, climate change, and sampling stage on the collembolan abundance and diversity. Our results suggest that Bt maize had no significant effect on the soil collembolan community. This is in line with our hypothesis (1) that Bt maize would have no significant effect on the soil collembolan community. Similarly, this is consistent with the findings of our previous study undertaken on a larger field (10 m × 15 m plot), in which we found that maize IE09S034 cultivation had no impact on the soil collembolan community 45 . Moreover, redundancy analysis also showed that the collembolan composition was not significantly correlated with the maize type. All these results indicated that Bt maize IE09S034 did not have adverse effects on soil Collembola. Although the results of previous studies have shown that the effect of Bt maize on soil Collembola varied with maize type and study sites, some studies have found that some species 46  Interestingly, we found that the climate changes simulated by OTC had a significant effect on the soil collembolan communities, and the effects varied by the sampling stage. We showed that the soil humidity was significantly lower and the soil temperature was significantly higher in the OTC than outside at all sampling stages. Temperature and humidity are the two most important climate variables, and changes in these can affect the community structure and ecological function of soil fauna by affecting their growth, the environmental conditions of their habitat, and by regulating their "bottom-up" food sources 52,53 . The abundance and other diversity indices of Collembola were significantly higher in the OTC than outside at the seedling stage, but significantly lower in the OTC than outside at the heading and mature stages. This may be related to the physiological characteristics of Collembola, the seasonal changes in temperature and humidity in northeast China, and the root biomass of maize. The exoskeleton of soil Collembola is thin and poorly sclerotized 54 ; thus, they are sensitive to soil temperature and humidity 28,29,55 . Our redundancy analysis also demonstrated that the collembolan composition was significantly correlated with the soil humidity and root biomass. In northeast China, the temperature was very low, usually < 0 °C before April. Under these conditions, the soil Collembola were basically dormant. Subsequently, as the temperature gradually rised, the soil Collembola began to revive and become active. Until the maize seedling stage (21th June, ~ 25 °C), the soil Collembola were in the recovery period. During this period, because the temperature inside the OTC was higher than that outside, soil Collembola recovered faster in the OTC than outside. Moreover, at this stage, the maize plant was still small, and the root biomass in the OTC was similar to that outside the OTC; therefore, the root exudate, which is the main food source of Collembola, was also not significantly different. Hence, the abundance and diversity of soil Collembola in OTC were higher than those outside OTC. Subsequently, from the seedling stage to the heading stage (4th August), the temperature continued to rise to ~ 35 °C, which is higher than the optimal living temperature for Collembola 56 , and the humidity in OTC (~ 25%) became significantly lower than the outside (~ 35%). This restricted the growth and development of Collembola 55,57 . Further, at this stage, the root biomass in the OTC was less than that outside. Therefore, the collembolan abundance and diversity in the OTC at the heading stage were significantly lower than those outside. At the mature stage (8th October), although the temperature had gradually decreased to the suitable temperature for Collembola (~ 18 °C), the humidity (~ 11%) and root biomass inside the OTC were still low. Consequently, along with the lower abundance and diversity of Collembola in the OTC at the heading stage, their reproduction was probably less than that outside; thus, the collembolan abundance and diversity in OTC were still significantly lower than those outside. Our results indicate that warming combined with drought influenced soil Collembola as follows: when the temperature is low, the soil collembolan abundance and diversity will be promoted by warming, and when the temperature is high and the humidity is low, soil collembolan abundance and diversity will be inhibited by warming combined with drought. Previous studies have also shown that the response of soil fauna to warming and drought depended on the ambient climate conditions 58 . For example, Harte et al. 59 found that warming increased the micro arthropod biomass in a cold, wet summer but reduced the micro arthropod biomass in the following warm, dry summer. Climate change resulting in frequent summer droughts would probably decrease the abundance and diversity of forest soil fauna 60 , while extreme winter warming events have been shown to favor large-bodied, litter-dwelling Collembola 61 .
However, no effects of the interactions of the maize type with the OTC were detected on these parameters. Some previous studies have found that warming and drought can alter the ingredients and quantity of plants 16,24 , which may, in turn, affect the underground ecosystem by changing the organic matter input. Here, we found that the root biomass of maize was significantly altered when the ambient conditions changed to reflect those anticipated with climate change. However, this alteration had no relationship with the maize type, and no effects of the interactions of the maize type with the OTC or the maize type with the OTC and the sampling stage on Collembola were found. This result conflicts with our hypothesis that climate changes simulated by OTC would alter the effects of Bt maize on soil Collembola, indicating that warming combined with drought did not alter the effect of Bt maize on soil Collembola.
In summary, we demonstrated that the effects of climate change on soil Collembola depended on the ambient climatic conditions due to the thin and poorly sclerotized exoskeleton of soil Collembola: when the temperature is low, Collembola are promoted by warming, and when the temperature is high and the humidity is low, Collembola are inhibited by warming combined with drought. Bt maize had no significant effect on the soil Collembola, and the effect could not be altered by the climate change simulated by OTC.  A distance of 20 cm between two adjacent maize plants in the same ridge, and a distance of 50 cm between two adjacent ridges were kept. Adjacent plots contained maize of different types. The same plot was used for planting the same maize type for the 3 years. No insecticides were applied during the study.
A humidity/temperature automatic recorder (Em50 Digital Data Logger, Decagon Devices Inc., Pullman, WA, USA) was placed in each plot to monitor the soil temperature and humidity of the 5 cm layer (recorded once per hour) in 2017. On each sampling day, at each sampling time (6 A.M., noon, and 6 P.M.), the soil temperature and humidity were measured using temperature and humidity recorders, and the mean of the three sampling times was used as the final data of the day.
Sample collection. Samples were collected at three crop stages in 2017: seedling (21th June), heading (4th August), and mature stages (8th October). Four sampling points within each plot were randomly selected for each sampling time, from where ~ 200 mL of soil was collected in a plastic bag using a soil auger (15 cm in diameter, 10 cm in height). Simultaneously, the maize plants near the sampling points were dug out and the roots were cut off. Following this, the soil and the roots were transported to the laboratory.
The 200 mL soil sample was placed on a 20-mesh screen over a Macfadyen extractor funnel 62 , and the soil Collembola were extracted for 7 days at room temperature (20-25 °C). The soil Collembola moved down through the funnel and dropped into a collection bottle below containing 95% ethanol. The number and species of the Collembola were counted using a microscope (Olympus SZ51, Japan) at a magnification of about 40×. The Collembola were identified to species based on previously published literary sources 63-66 using a microscope (Nikon 80i, Japan) at a magnification of about 200×. The roots were washed and dried in an oven at 55 °C for 48 h. The root biomass was then measured and the mean biomass of the four roots from one plot was used as the final data of the plot.

Statistical analysis. Dominance (D) of soil Collembola is the proportion of individuals belonged to
the ith species in the samples of the same treatment. Species with D > 10% were considered as numerically dominant 67,68 . The collembolan abundance and diversity parameters (species richness, Shannon-Wiener index, and Pielou's evenness index) were calculated using Data Processing System software (version 2005, China). The Shannon-Wiener index (H´) and Pielou's evenness index (J) were calculated as follows: where P i is the proportion of individuals belonged to the ith species in one sample.
where S is species number of the collected Collembola in one sample.
Three-way ANOVA was performed using SPSS (version 23, IBM, USA) to evaluate the effects of the maize type, OTC treatment, and sampling stage on the soil collembolan abundance and diversity (species richness, Shannon-Wiener index, and Pielou's evenness index). The collembolan abundance was log 10 (N + 1) transformed before the analyses to ensure normality and equal variances and Levene's test of equality of error variances and asymmetry were used to analyze the homogeneity of variance. Tukey's HSD post hoc test (α = 0.05) was used when significant differences between treatment means occurred. Moreover, a one-way ANOVA was performed to evaluate the effects of OTC on the soil collembolan abundance and diversity, and certain important environmental factors (soil humidity, soil temperature, and maize root biomass) for the same maize type at each sampling stage using SPSS.